Entanglement Sudden Death in Band Gaps 
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Using the pseudomode method, we evaluate exactly time-dependent entanglement for two independent qubits, 
each coupled to a non-Markovian structured environment. Our results suggest a possible way to control entan- 
glement sudden death by modifying the qubit-pseudomode detuning and the spectrum of the reservoirs. Par- 
ticularly, in environments structured by a model of a density-of-states gap which has two poles, entanglement 
trapping and prevention of entanglement sudden death occur in the weak-coupling regime. 
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I. INTRODUCTION 



Realistic quantum systems are affected by decoherence and 
entanglement losses because of the unavoidable interaction 
with their environments [1]. For example, in Markovian 
environments, in spite of an exponential decay of the single- 
qubit coherence, the entanglement between two qubits may 
completely disappear at a finite time [2,3]. This phenomenon, 
known as entanglement sudden death and experimentally 
proven to occur [4,5], limits the time when entanglement 
can be exploited for practical purposes. So the issue of 
how to avoid or control entanglement sudden death-type 
decorrelation in a realistic physical system is especially im- 
portant. Up to now, several ways were proposed for keeping 
the atomic entanglement for a long time by suppressing 
spontaneous emission. One way widely applied is to place 
the qubits in a structured environment, say, microcavity [6, 7] 
or in the photonic band gap of photonic crystals [8]. Other 
methods considered use dynamic manipulation such as mode 
modulation [9] or the quantum Zeno effect [10], which can 
be regarded as extensions of the so-called bang-bang method 
[11]. 

In this paper we focus on the entanglement dynamics 
of two qubits, each coupled to a non-Markovian structured 
environment. We continue the investigation of physical 
effects that may control the occurrence time of entanglement 
sudden death, and find that the speed of disentanglement 
is closely related to qubit-pseudomode detuning and the 
spectrum of the reservoirs. In particular, we mainly consider 
two qubits respectively coupled to a bath of oscillators with a 
density of states described by a frequency dependent function 
D(ui). In some previous works [12, 13, 14], the spectrum 
of the structured reservoir is a Lorentzian, which only give 
a single pseudomode to replace the effect of the structured 
reservoir. Then the short-time behavior of the qubit system 
leads to exponential decay [15]. However, if there are two or 
more poles close to the real cu axis, it is not clear that a simple 
picture of exponential decay will apply because the poles 
can interfere with each other. So we focus on the structured 



reservoir model of a band gap in which both Lorentzians are 
centered at the same frequency, and the second is given a 
negative weighting in our paper, and make comparison with 
the model of the function D(uS) has one Lorentzian. 



II. THEORETICAL MODEL AND EXACT DYNAMICS OF 
TWO QUBITS 

Consider a model consisting of two qubits A and B, each 
interacting with a common zero-temperature bosonic reser- 
voir, denoted a and b, respectively. We assume that each 
qubit-reservoir system is isolated and the reservoirs are ini- 
tially in the vacuum state while two qubits are initially in an 
entangled state. Since each qubit evolves independently, we 
can learn how to characterize the evolution of the overall sys- 
tem from the qubit-reservoir dynamics. The interaction be- 
tween a qubit and an iV-mode reservoir is described through 
the Hamiltonian (under the rotating-wave approximation and 
setting h — 1) 
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where bl, bk are the creation and annihilation operators of 
quanta of the reservoir (a or b), a J + = |1j)(0j|, <3~1 = |0j)(l 3 -| 
and u>j are the inversion operators and transition frequency of 
the j-th qubit Q=A, B and here loa = o->b = )', w fe an d 9k 
are the frequency of the mode k of the reservoir and its cou- 
pling strength with the j-th qubit. 

Let us consider the case when the j-th qubit is initially in the 
excited state and its corresponding reservoir is in the vacuum 
state. Using the pseudomode method, we focus on an ideal- 
ized model of a band gap (or photon density of states gap) in 
which both Lorentzians are centered at the same frequency, 
and the second is given a negative weighting, so that [16] 
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where now the weights of the two Lorentzians are such that 
Wi — W2 = 1 and T 2 < Ti ensure positivity of D. uj c is 
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the center of the spectrum, and Ti, T 2 are the full width at 
half maximum of two Lorentzians, respectively. The effect 
of the Lorentzian with negative weight is to introduce a dip 
into the density of states function D(u>) where the coupling 
of the qubit will be inhibited [17, 18]. The two poles in D are 
located at uj c — iT i/2 and lj c — iT 2 /2, and there is a change 
in sign of the residues of D between these poles. The exact 
pseudomode master equation associated with the density of 
states gap in Eq.© is given by 

dp r' 

— = - i[H J ,p] - ylajaip- 2a 1 pa{ + pa^} 

r' 

- -yt a 2 a2 ^ - 2a 2 pa\ + pa\a 2 ], (3) 

where 

Hi = ujQa-^a^ + u) c a\ai + uj c a\a 2 + ilo[a 2 o"i + «20+] 
+ V(a\a 2 + axal), (4) 

where p is the density operator for the j-th qubit and the pseu- 
domodes of its corresponding reservoir; and a\ and a 2 are the 
annihilation operators of the two pseudomodes decaying with 
deacy rates Ti = W X T 2 - W 2 T 1 and T' 2 = WiTi - W 2 T 2 
respectively. The two pseudomodes are coupled and V = 
y/WiW 2 (Ti — T 2 ) /2 is the strength of the coupling. The qubit 
interacts coherently with the second pseudomode (the strength 
of the coupling Oo), which is in turn coupled to the first one. 
Both pseudomodes are leaking into independent Markovian 
reservoirs. The set of ordinary differential equations associ- 
ated to the master equation ([3]) is 

dc\ 

1— = loqCx + Q a b 2 , 
at 

■ dbl >h LIM 

1— = z x b\ + Vb 2 , 
at 

i-rr = z 2 h + Vbi+Q, ci, (5) 
at 

where c±, b\, and b 2 are the complex amplitudes for the 
states with one excitation in the qubit, one excitation in the 
first pseudomode, and one excitation in the second pseu- 
domode, respectively. The positions of the true poles are 
z[ = u) c — iT^/2 and z' 2 — oj c — «r 2 /2. Due to the ini- 
tial state of the j-th qubit and its corresponding reservoir is 
® |0> (with the state |0) = Uk=i \°k)), then c x (0) = 1, 
and 6i(0) = 62(0) = 0, so we can acquire the exact solutions 
(ci, b\ and b 2 ) of the differential equations (0 easily through 
a computer program. 

According to the differential equations (|5), for an initial 
state of the form 

|*(0)) = (a|00) i4 B + /3|ll)>ifl)® |0>„|0) 6 , (6) 

where a, (3 are real, then the time evolution of the total system 
is (using the above pseudomode method) 

|*(t)) = /3(c 1 |l) A |0 1 2 ) a + 6 1 |0) A |li0 2 > a + 6 2 |0>A|0 1 l 2 ) a ) 
x (ci|l) B |0i0 2 ) b + &i|0)b|1i0 2 > 6 + fe 2 |0) s |0il 2 ) b ) 
+ a|00) A B|0i0 2 ) a |0i0 2 )6, (7) 



with |0i0 2 ) a /f„ |li0 2 ) a /b and |0il 2 ) a /b mean the structured 
reservoir a or b states with no excitation in two pseudomodes, 
only one excitation in the first pseudomode, and one excitation 
in the second pseudomode, respectively. Then we can deter- 
mine the two qubits dynamics. In particular, in the standard 
two qubits basis C = {|00)ab, |01)ab, \W)ab, \H)ab}, the 
reduced density matrix of the two qubits at time t result as 

PA B (t) = (a 2 +/3 2 (l- |ci| 2 ) 2 )|00)ab<00| 

+ (/3 2 | C1 | 2 (1-| C1 | 2 ))|01)^<01| 

+ (/3 2 | Cl | 2 (l-| Cl | 2 ))|10)^<10| 

+ /3 2 | C1 | 4 |11) AB <11| 

+ a/?c 2 |00) AB (ll| 

+ a(3{c\y |11) AB (00|. (8) 

We use the concurrence C [19], which attains its maximum 
value 1 for maximally entangled states and vanishes for sepa- 
rable states, to analyze the two-qubit entanglement dynamics. 
For pAB{t), its concurrence can be derived from [19], as 

C{t) = 2max{0,a/3| Cl | 2 - /3 2 | Cl | 2 (l - | Cl | 2 )}. (9) 

III. NUMERICAL RESULTS AND DISCUSSIONS 

Similar to the result in Ref. [20], it is easy to find that the 
two-qubit entanglement can occur sudden death for a < (3. 
In this non-Markovian system-reservoir coupling model, the 
concurrence of the two qubits will vanish forever at a finite 
interval in weak-coupling regime, but for the strong-coupling 
regime, there will be entanglement revival after entanglement 
sudden death. In the following, we mainly investigate that 
the speed of the occurrence of entanglement sudden death is 
related to the spectrum of the reservoirs (T 2 ) and the qubit- 
pseudomode detuning (A = ui c — u>q). 

First we study the relation between the entanglement sud- 
den death and the spectrum of the reservoirs. Fig.l shows the 
entanglement dynamics of two qubits in the non-Markovian 
weak-coupling regime with Ti/2 = 10£!o- We compare the 
entanglement dynamics of the two qubits for three different 
values of the width of the second Lorentzian spectral function, 
namely, T 2 /2 = ilo, 2Hq, 9f2o- As in Fig. 1(a), the qubits are 
on resonance with the center of the spectrum, A = 0. The 
concurrence decreases monotonically down to zero in a finite 
time. It is interesting to find that the speed of occurrence of 
entanglement sudden death can increase with T 2 increasing. 
Similar behavior happens when two qubits are near resonance 
with the center of the spectrum. However, When the qubits 
are far off -resonant with the center of the spectrum, A>Oo, 
as shown in Fig. 1(6), where we choose A = lOf^o, the speed 
of occurrence of entanglement sudden death decreases as T 2 
increases. We can give an intuitive explanation for these re- 
sults. As we can see in Fig. 2, the density of the spectrum 
D(lu) increases monotonically as T 2 increases for A = 0, 
while the density of the spectrum D(u) decreases monotoni- 
cally as r 2 increases when A = 10f2o. It is proved that the 
entanglement sudden death is determined by the modes of the 
spectrum which are on resonance with the qubits: the speed 
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FIG. 1 : Time evolution of the concurrence of two qubits as a function 
of the dimensionless quantity flot in non-Markovian weak-coupling 
regime, with (a) A = and (b) A = lOfio- For the cases of (i) 
T2 /2 = Oo (solid dark curve), (ii) F2 /2 = 2Qo (dashed red curve), 
(iii) T2/2 = 9f2o (dotted blue curve). The parameters used are: 
Ti/2 = lOfio, Wi = 1.1, W 2 = 0.1 and a = §,/3 = ^. 



FIG. 2: The density of the spectrum D(u>) as a function of the dimen- 
sionless quantity 1^/2 in non-Markovian weak-coupling regime, 
with (a) A = and (b) A = lOfio- The parameters used are: 
Ti/2 = lOfio, Wi = 1.1, W 2 = 0.1. 



of the occurrence of entanglement sudden death decreases (in- 
creases) as the density of these modes decreases (increases). 

If the structured reservoir only contains a Lorentzian 
(T2 = 0), from the Ref. [12], we note that the speed of disen- 
tanglement decreases as the width of the Lorentzian spectral 
(Ti is replaced by A in Ref. [12]) increases on/near the reso- 
nant couplings, and the speed of disentanglement increases as 
A increases for large deduning coupling. However, in our pa- 
per the density of the spectrum D(w) is composed by two 
Lorentzians, and the second is given a negative weighting. 
Then the counter results can be acquired when we consider 
the width of the second Lorentzian spectral T2 influence the 
speed of the occurrence of entanglement sudden death. That 
is to say, Ti and T2 have the opposite effects on the speed of 
the occurrence of entanglement sudden death. So we can keep 
two-qubit entanglement for a long time trough choosing suit- 
able spectrum of the reservoirs. 

As we know, the effect of the Lorentzian with negative 
weight is to introduce a dip into the density of states func- 
tion D(ui). For a perfect gap, where D(lo c ) = 0, we 
would also have W1/T1 — W2/T2, then the single-qubit 
excited-state population trapping can occurs [15]. According 
to Ref. [21], we can find that there is a direct link between 
the time-dependent entanglement and single-qubit excited- 
state population for independent qubits, each coupled to a 
zero-temperature bosonic environment. So there also will ap- 
pear two-qubit entanglement trapping if the qubits are reso- 



nant with the gap in our paper. But it is surprising to see 
that the entanglement trapping only can occur in the weak- 
coupling regime, and entanglement sudden death still ap- 
pears quickly in the intermediate-coupling regime and strong- 
coupling regime. For the intermediate-coupling regime, the 
two-qubit entanglement can vanish forever. Due to the strong 
non-Markovian effects, the entanglement of the two qubits 
can arise the phenomenon of entanglement sudden death and 
revival in the strong-coupling regime (as shown in Fig. 3). 

Then we fix the spectrum of the reservoirs to study the 
relation between the entanglement sudden death and qubit- 
pseudomode detuning A. In Fig. 4 we show the entangle- 
ment dynamics for two qubits in a structured reservoir, which 
only has one Lorentzian spectral. Fig. 5 shows the entangle- 
ment dynamics when the qubits interact with two indepen- 
dent reservoirs structured by a model of a density-of-states 
gap which has two Lorentzians. A comparison between Fig. 4 
and Fig. 5 reveals that, both for the weak-coupling regime, 
the speed of the occurrence of entanglement sudden death 
only can decrease as |A| increases in the one Lorentzian 
model. However, there exist a critical value |A C | (with A;? = 

r ^J^^J^-^J^y - ) m me band gap model, then the speed of 
occurrence of entanglement sudden death will increase as | A| 
increases when |A| < |A C |, and entanglement sudden death 
can occur more slowly with |A| increasing at |A| > |A C |. 
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FIG. 3: Time evolution of the concurrence of two qubits as 
a function of the dimensionless quantity iloi, for the cases of 
(i) weak-coupling regime ri/ 2 = Hf2o,r2/2 = flo (solid dark 
curve), (ii)intermediate-coupling regime Ti/2 = 1.1X20)^/2 = 
O.lfio (dashed red curve), (iii)strong-coupling regime Ti/2 = 
0.11f2o,r2/2 = O.Olfio (dotted blue curve). The parameters used 
are: Wi = 1.1, W 2 = 0.1 and a = \,p= ^. 




FIG. 6: The density of the spectrum D(u) as a function of the dimen- 
sionless quantity A in non-Markovian weak-coupling regime, for the 
one Lorentzian model Ti/2 = 10f2o, Wi = 1 (solid dark curve) 
and the band gap model Ti/2 = lOfto, T 2 /2 = f2 , Wi = 1.1 and 
Wi = 0.1 (dashed red curve). 



A physical interpretation of the result is that the density of 
the spectrum D(uj) decreases monotonically as |A| increases 
in the one Lorentzian model, while in the band gap model, 
D(ui) increases as |A| increases when |A| < |A C | and de- 
creases as |A| increases at |A| > |A C |. The relation between 
D(oj) and A is clearly shown in Fig. 6, for example, when 
Ti/2 = lOfio, T 2 /2 = fl , we can calculate A c = ±3.53fi - 



IV. CONCLUSION 



FIG. 4: Time evolution of the concurrence of two qubits as a function 
of the dimensionless quantity £lot and | A| in reservoirs structured by 
a model of one Lorentzian spectral density. The parameters used are: 
Wi = 1, Ti/2 = lO^o and a = \,(3 = ^. 




FIG. 5: Time evolution of the concurrence of two qubits as a function 
of the dimensionless quantity Q,ot and |A| in reservoirs structured 
by a model of a density-of-states gap which has two Lorentzians. 
The parameters used are: Wi = 1.1, W 2 = 0.1, Ti/2 = 10fi , 
T 2 /2 = fi , and a = \,fi=^. 



In summary, we have presented a non-Markovian model de- 
scribing the exact entanglement dynamics of two qubits, each 
interacted with a structured reservoir. We have brought to light 
new relations between the speed of the occurrence of entan- 
glement sudden death and the spectrum of the reservoirs/the 
qubit-pseudomode detuning for qubits prepared in entangled 
state. We firstly find that the speed of occurrence of entan- 
glement sudden death is a increasing (decreasing) function 
of the width of the second Lorentzian spectral T2 when the 
qubits are on/near resonance (large detuning) with the pseu- 
domodes of the reservoirs. Due to the density of the spectrum 
D(uj) is composed by two Lorentzians in our paper, and the 
second is given a negative weighting, it is interesting to find 
that Ti and T2 have the opposite effects on the speed of the 
occurrence of entanglement sudden death. Then for a per- 
fect gap, where D(oj c ) — 0, we would also have Wi/Ti = 
W2/T2, the qubits' entanglement trapping and prevention of 
entanglement sudden death can occur in the weak-coupling 
regime, but entanglement sudden death still appears quickly in 
the intermediate-coupling regime and strong-coupling regime. 
Next, a comparison between the one Lorentzian model and the 
band gap model, we find the speed of the occurrence of entan- 
glement sudden death only can decrease as |A| increases in 
the one Lorentzian model and a critical value |A C | exist in 
the band gap model. If |A| < |A C |, the speed of occurrence 
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of entanglement sudden death will increase as |A| increases, 
and when |A| > | A c |, entanglement sudden death can occur 
more slowly with |A| increasing. The results here obtained 
evidence the entanglement can be preserved or controlled by 
modifying the spectrum of the environment and highlight the 
potential of reservoir engineering for controlling and manipu- 
lating the dynamics of quantum systems. 

Our results would apply to cavity QED experiments with 
trapping ions, and to circuit QED experiments. Entanglement 
between two remotely located trapped atomic ions has been 
recently demonstrated [22] and multiparticle-entangled states 
can be generated and fully characterized via state tomogra- 
phy [23]. Moreover, field coupling and coherent quantum 
state storage between two Josephson phase qubits has been 



achieved through a microwave cavity on chip [24]. Due to the 
possibilities for realizing strong coupling conditions between 
atoms and a high finesse cavity [25], a deep understanding of 
the non-Markovian dynamics is now indispensable. 
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